SEARCH AND TEST ALGORITHMS FOR 
TRIPLE PRODUCT PROPERTY TRIPLES 



IVO HEDTKE AND SANDEEP MURTHY 



Abstract. In 2003 Cohn and Umans introduced a group-theoretic approach to fast matrix 
multiplication. This involves finding large subsets of a group G satisfying the Triple Product 
Property (TPP) as a means to bound the exponent uj of matrix multiplication. 

We present two new characterizations of the TPP, which are useful for theoretical consider- 
ations and for TPP test algorithms. With this we describe all known TPP tests and implement 
them in GAP algorithms. We also compare their runtime. Furthermore we show that the search 
for subgroup TPP triples of nontrivial size in a nonabelian group can be restricted to the set 
of all nonnormal subgroups of that group. Finally we describe brute-force search algorithms 
for maximal subgroup and subset TPP triples. In addition we present the results of the sub- 
set brute-force search for all groups of order less than 25 and selected results of the subgroup 
brute-force search for 2-groups, SL„ Fq and PSL2Fq. 

1. Introduction 

1.1. A Short History of Fast Matrix Multiplication. The naive algorithm for matrix mul- 
tiphcation is an Olrr') algorithm. From Volker Strassen (see [Strassen 1969) ) we know that 
there is an 0{1l?^^^) algorithm for this problem. Winograd optimized Strassen's algorithm. 
While the Strassen- WiNOGRAD algorithm is the variant that is always implemented (for ex- 
ample in the famous GEMMW package), there are faster ones (in theory) that are impractical 
to implement. The fastest known algorithm runs in 0{'}l?^^^) time (see [Coppersmith et al. 1987| 
from Coppersmith and Winograd). Most researchers believe that an optimal algorithm with 
0{n^) runtime exists, but since 1987 no further progress was made in finding one. Because modern 
microprocessor architectures have complex memory hierarchies and increasing parallelism, perfor- 
mance has become a complex tradeoff, not just a simple matter of counting flops. Algorithms 
which make use of this technology were described in |D 'Alberto et al. 2009) . Another well known 
method is Tiling: The normal algorithm can be speeded up by a factor of two by using a six loop 
implementation that blocks submatrices so that the data passes through the LI Cache only once. 

1.2. The Exponent of the Matrix Multiplication. Let M{n) denote the number of field 
operations in characteristic required to multiply two {n x n) matrices. The exponent uj of the 
matrix multiplication is defined as 

LO := inf{r G M : M{n) = 0{rf)}. 

Details about the complexity of matrix multiplication and the exponent uj can be found in 
Biirgisser et. al 1997] . The result from Coppersmith and Winograd says that uj < 2.38. 

1.3. The Group-theoretic Approach from Cohn and Umans. (The reader can find the nec- 
essary background on group- and representation theory in jAlperin et al. 19"9T] , [Jame s et al. 2011] 
and [Neumann 1994] . More details about the framework can be found in [Cohn et al. 2003) and 
[Cohn et al. 2005] .) 

In 2003 Cohn and Umans introduced in [Cohn et al. 2003] a group-theoretic approach to fast 
matrix multiplication. The main idea is to embed the matrix multiplication over a ring R into the 
group ring RG, where G is a (finite) group. A group G admits such an embedding, if there are 
subsets S, T and U which fulfill the so-called Triple Product Property. 
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Definition 1.1 (right quotient). Let G be a group and 7^ X C G be a nonempty subset of G. 
The right quotient Q{X) of X is defined by Q{X) := {xy^^ : x,y ^ X}. 

Note that Q{S) = S holds, if S* is a subgroup of G. 

Definition 1.2 (Triple Product Property). We say that the nonempty subsets S, T and U of 
a group G fulfill the Triple Product Property (TPP) if for s e Q(5), t e Q{T) and u G Q([/), 
stu — 1 holds iffs = t = u = l. 

With {n,p,m) we denote the problem 

{n,p,m): C"^p x C^^" ^ C"''", (^, B) AS, 

of multiplying an (n x p) with a (p x to) matrix over C. We say that a group G realizes (si, S2, 53) 
if there are subsets S'i C G of sizes \Si\ = Si, which fulfill the TPP. In this case we call (^i, 52, S3) 
a TPP triple of G. Let us now focus on the embedding of the matrix multiplication into CG. Let 
G realize {n,p,m) through the subsets |5| — n, \T\ = p and \U\ = to. Let A be an (n x p) and B 
be a (p X to) matrix. We index the entries of A and B with the elements of S, T and U instead 
of numbers. Now we have 

{AB)s^u^J2t^T^'-'^''^- 

COHN and Umans showed that this is the same as the coefficient of s^^u in the product 

So we can read off the matrix product from the group ring product by looking at the coefficients 
of s~^u with s G S and u G U. 

Definition 1.3 (TPP (subgroup) capacity). We define the TPP capacity /3(G) of a nontrivial 
group G as /?(G) := maxjnpm : G realizes {n,p,m)} and the TPP subgroup capacity of G as 
/3g(G) :— maxjnpTO : G realizes {n,p,m) through subgroups}. 

Note that /3(G) > |G|, because every group G realizes (|G|, 1, 1) through the TPP triple (G, 1, 1). 

Definition 1.4 (r-character capacity). Let G be a group with the character degrees {di}. We 
define the r-character capacity of G as Dj,{G) :— 

We can now use (3 and Dr to get new bounds for uj: 

Tlieorem 1.5. jCohn et al.3003l Thm. 4.1] If G ^ 1 is a finite group, then l3{Gy/^ < D^{G). 

Note that this leads to a nontrivial upper bound for w, iff /3(G) > Dy,{G). Therefore, for a 
fixed group G we search for TPP triples (S, T, U) which maximize |5| • |r| • \U\, for example with 
a brute-force computer search. Note that the maximal \S\ ■ \T\ ■ \U\ equals to /3(G). 

With their group theoretic framework Cohn and Umans were able to construct an algorithm 
for matrix multiplication with running time ©(n^'^^), see |Cohn et al. 2005"! Sec. 6.3]. 

From COHN and Umans we know f |Cohn et al. 20031 Lem. 3.1.]), that ^(G) = |G| if G is an 
abelian group. Therefore we only focus on nonabelian groups. We also only focus on finite groups. 

1.4. Tiie Aim of tiiis Work. In this article we present two new characterizations of the TPP, 
which are useful for theoretical considerations and for TPP test algorithms. With this we describe 
all known TPP tests and implement them in GAP. We also compare their runtime. Furthermore 
we show that the search for subgroup TPP triples of nontrivial size in a nonabelian group can 
be restricted to the set of all nonnormal subgroups of that group. Finally we describe brute- 
force search algorithms for maximal subgroup and subset TPP triples. In addition we present the 
results of the subset brute-force search for all groups of order less than 25 and selected results of 
the subgroup brute- force search for 2-groups, SL„ and PSL2 F^. 
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2. Basic Concepts 

In this section we present facts that we use in our search algorithms. Details about the subset and 
subgroup search can be found in section [6] Depending on the context, the symbol 1 will denote 
either the number 1, the group identity Iq, or the trivial subgroup {1g}- 

Lemma 2.1. [Hedtke 20lTl Lem. 2 and 3 and Thm. 4] 

(*) Let $ X C G be a nonempty subset of a group G and g € G. Then 1 G Q{X) and 
geQ{X)^g~^ eQ{X). 
{**) If S, T and U fulfill the TPP then Q{X) n Q{Y) = 1 holds for allX e {S, T, U}. 
(***) // {S, T, U) is a TPP triple with l€ SCiTnU, then S f] T ^ T f] U ^ S H U ^ 1. 

Lemma 2.2. [Hedtke 20111 Cor. 6] If{S,T,U) is a TPP triple ofG, then\S\ + \T\ + \U\ < |G|+2. 

First we note, that we are only interested in matrix-matrix multiplication, that means, we 
consider only TPP triples {S,T,U) with 15*1, |T|, \U\ > 1. Furthermore we have: Assume that in 
a TPP triple of G one of S, T or U is G itself. Then it follows from the lemma above, that the 
other two sets of the triple have size 1. Because we omit this case, we only use: 

Observation 2.3. It is sufficient to search TPP triples with \S\, \T\, \U\ e {2,...,\G\~ 1}. 

From Corn and Umans (see |Cohn et al. 20031 Lem. 2.1]) we know the following symmetry 
property: 

Lemma 2.4. If G realizes {n,p,m), then it does so for every permutation of n, p and m. 

We conclude the following: 
Observation 2.5. It is sufficient to search TPP triples with \S\ > \T\ > \U\. 

We know, that /3(G) > |G|, therefore we are only interested in 

Observation 2.6. It is sufficient to search TPP triples with \S\ ■ \T\ ■ \U\ > \G\. 

Now we can combine the last two observations: We know that \S\ ■ \T\ ■ \U\ < \S\'^. Therefore 
we have 

Observation 2.7. It is sufficient to search TPP triples with \S\ > y/\G\ —: i{G), where \x~\ 
denotes the smallest integer not smaller than x. 

Assume the case, where we try to find triples that yield a nontrivial upper bound for w instead 
of only finding nontrivial TPP triples. In this case we change the |G| in the last two observations 
to D^[G). 

Definition 2.8 (basic TPP triple). According to Neumann we call a TPP triple {S,T,U) that 
fulfills leSTlTnt/a basic TPP triple. 

From Neumann we know the following facts that reduce the search space enormously: 

Lemma 2.9. (Neumann 20111 Obs. 2.1] If {S,T,U) is a TPP triple of G, then {dSa,dTb,dUc) 
is a TPP triple for all a, b,c,d€ G, too. 

Note that any TPP triple can be translated to a basic TPP triple by lemma [279] Together with 
(***) it follows, that: 

Observation 2.10. It is sufficient to search TPP triples with S (IT ^ T nU ^ S nU = 1. 

Lemma 2.11. [Neumann 20111 Obs. 3.1] If{S,T,U) is a TPP triple, then |S'|(|T| + |t/|-l) < |G|, 
md^l + \U\-1)<G and \U\i\S\ + \T\ - 1) < |G|. 

Observation 2.12. It is sufficient to search TPP triples with \S\i\T\ + \U\ - 1) < |G|. 

Example 2.13. Let us focus on a brute-force search for TPP triples of subgroups in the Mathieu 
group Mil. We pick out this example because it is a big group, note that |A/ii| = 7920, with 
many (exactly 8651) subgroups. Possible orders for S, T or U are (computed with GAP): 

{1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 16, 18, 20, 24, 36, 48, 55, 60, 72, 120, 144, 360, 660, 720, 7920} =: V. 
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Note that V has 26 elements. If we use observation 2.3 we define V := V \ {1,7920}, so we 
have \r\ = 24. A naive idea would be {\Sl\T\,\U\) e =: C. Because \C\ = 13824 we 
would have a large search space for a brute-force search. Now we use observation |2.5| and achieve 
|C| ~ 2600. If we only search for subgroup TPP triples that yield a nontrivial bound for uj, we 
use observation |2.6| with D^{G) = 355208 (again computed with GAP) instead of \G\. With this 



it follows that |C| = 404. Finally we use Neumanns inequality from observation 2.12 which 
yields to |C| — 0. This example shows that the observations above reduce the search space for a 
brute-force search enormously. In this case it is not even necessary to start a search. 

3. Two New Characterizations of the Triple Product Property 
In this section we present two new characterizations of the TPP. They are useful for theoreti- 



cal considerations (in particular theorem 3.1) and for TPP test algorithms which we discuss in 
section [5l 

Theorem 3.1. Three subsets of G form a basic TPP triple {S,T,U) iff 

(i)le5nTnC/, (ii) Q(T) n Q([/) = l and (iu) Q{S) n Q{T)Q{U) = 1. 

This is not a limitation, because we only need to search for basic TPP triples. 

Proof. First assume that {S, T, U) is a basic TPP triple, (i) follows directly from the definition and 
(ii) from (**). Furthermore (*) implies that 1 S Q{S)r)Q{T)Q{U). Now assume there is a common 
element 1 ^ x e Q{S) n Q{T)Q{U). Then x = s = tu, for some s € Q{S) and tu e Q{T)Q{U). 
This means 1 — x^^x = s^^tu, but the TPP for {S, T, U) implies that 1 = s^^ = t = u and 
therefore x = 1, a contradiction. 

Now assume that the equations (i)-(iii) hold for the subsets S, T and U. Consider the triple 
quotient product stu for arbitrary elements s e Q{S), t G QiT) and u G Q{U). Then stu = 1 is 
equivalent to s"^ = tu. Now (iii) implies that = tu = 1 and (*) together with (ii) imply that 
t = u = 1 and so (5, T, U) is a basic TPP triple. □ 

Definition 3.2 (subtransversal, support). Let C be a finite nonempty set and C — {Ci, . . . , Ck} 
a partition of it. A set X Q C is called a subtransversal for C with support swppq{X) = T C C if 
for all e C 

fl QeT, 

I otherwise. 



IX n a 



It then follows that \X\ = |T|. In the special case when the collection C is the set of left (or right) 
cosets of a subgroup S* of a group G, then any subtransversal T for G/S (or S\G) will simply be 
called a subtransversal for S in G. 

Theorem 3.3. Let G be a group, S a subgroup of G, and T, U subsets of G. 

(1) // {S, T, U) is a basic TPP triple of G then T and U are subtransversals for S in G such 
that 

(1) snppg^ciT) n supps\GiU) = {S}. 

(2) IfT and U are also subgroups of G, and T and U are subtransversals for S in G satisfying 
g then {S,T,U) is a TPP triple ofG. 

Proof, (i) We have S r\T ^ S C\U ^ 1, hy (***). Let 5 7^ 5r e \ G be a nontrivial coset of S. 
Assume distinct elements 1 7^ t, t' G Tn Sr, where t = sr and t' = s'r for distinct s, s' e S. Then 
1 7^ t't^^ = s's~^ G S f] Q{T), and that contradicts the TPP requirement (**). So for any coset 
Sr E S \ G , we have |T n Sr\ < 1. The same holds for U. Therefore T and U are subtransversals 
for S \G. Now assume distinct elements l^t — srETDSr and 1 ^ u = s'r E U D Sr ior a 
nontrivial coset Sr ^ S. Then 1 ^ tu"! = s(s')"^ G Q(5) n Q{T)Q{U) = SC^ Q{T)Q{U) and that 
contradicts theorem [3J| This shows that T and U fulfill 



(ii) Assume T and U are subtransversals of 5'\G and fulfill Q. Since 5, T and U are subgroups 
we have 1 E S DT DU and 1 E S D TU. Because the intersection of the supports of T and U is 
{S}, TCiU = 1 holds. Now assume that there is an 1 7^ a; G 5 n TU. So there aie s E S, t eT 
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and u with 1 ^ s = tu. Therefore we have tu e S, which is equivalent to t G Su ^ . But then 



would be i = ?i = 1, a contradiction. So [S, T, U) is a TPP triple by theorem 3.1 □ 

Observation 3.4. Three subgroups S , T, U of a group G form a TPP triple iff T and U are 
subtransversals for S in G, satisfying ([ij . 

Theorem 3.5. Let G be a group. If {S, T, U) is a TPP triple of subgroups where at least one of 
S , T or U is normal in G, then \S\ ■ \T\ ■ \U\ < \G\. 

Proof. Without loss of generality assume that S is nontrivial, proper and normal in G. Because S 



is normal in G, we have S\G = G/S. We will work with G/S. From theorem 3.3 we know, that T 
and U are subtransversals for G/S that fulfill Let Si := 5, Sj, . . . , Sy be the v = [G : S] > 2 
elements of G/S. We define T suppQ/g{T) G/S and U in the same way. Then T and U are 
subgroups of G/S with |T| = \T\, \14\ — \U\ and T r\U = 1. To see this, first consider T. Because 
T and U are subtransversals for G/S that fulfill ([T]), we have Iq/s — S & T. Now let rS,r'S € T 
be cosets of S. So there are t,f e T such that T HrS ^ {t ^ rs} and T n r' S ^ {f = r's'}, 
where s, s' G S. They have the product rsr's' = tt' € T n rr' S and so \T n rr' S\ > 1. Since T 
is a subtransversal for G/S it follows that |T n rr'S\ = 1 and so rr'S G T. The inverse of any 
given rS € T is r~^S £ G/S and if t G T is such that THrS ^ {t} then t = rs for some s G S", 
and = s"^r^^ G S'r^^ — r~^S (because S is normal), and so t~-^ G T n r^^S, and we can 
conclude that r^^S G T as well. The same holds for U. Because ([l]) we know that T fMA = 1 



and so {TM, 1) is a TPP triple of G/S by theorem 3.1 From Neumanns inequality (see lemma 



2.9) it follows that |r|(|W| + 1 - 1) = |r| • \U\ = \T\ ■ \U\ < \G/S\ = \G\/\Sl which we wanted to 



show. □ 

We can use the result above to create an additional filter for the search space of subgroup TPP 
triples: 

Observation 3.6. // we are only interested in subgroup TPP triples of nontrivial size it is suffi- 
cient to search for S, T and U that are all nonnormal. 

4. On the TPP Capacity of Nonabelian Groups 

In this section we present some facts about the TPP capacity of nonabelian groups. The following 
two results state nontrivial lower bounds for /3(G). The ideas behind the proofs are those of the 
authors, but the proofs as formulated here are those of Neumann. 

Lemma 4.1. If G is a nonabelian group with a nonnormal subgroup S of index [G : S]— 3, then 

p{G)>m- 

Proof. Let S\G be the right coset space of G of size [G : S] = 3. There is a natural honiomorphism 
(j): G 5*3, defined hy g i-^ iTg for all g € G, where iTg describes a permutation action (g, Sr) 
Srg — Srg G 5 \ G for all iSr G 5' \ G, of g on S* \ G. The homomorphism has a kernel 
K := ker 0<G, which is CoreG'(S') (the largest normal subgroup of G contained in 5*). The quotient 
group G/K = img0 < S3 is isomorphic to a transitive subgroup in S3. The nontrivial transitive 
subgroups of S3 are A3 and 5*3 itself. Since S is nonnormal, K < S is proper and so \G/K\ > 3. 
Thus G/K = 5*3. From [Cohn et al. 20031 Lem. 2.2] we know that /3(G) > (i{G/K)P{K). It 
follows (see table[l[, that /3(G) > 8/3(/r) > 8|A'| = 8\G\/\G/K\ = ||G| = ||G|. □ 

Lemma 4.2. If G is a nonabelian group with a self-normalising subgroup S of index [G : S] = 4, 
then /3(G) > ||G|. 

Proof. The subgroup S has a normal core K :— GoiedS) < G which is proper K < S, such 
that the quotient group G/K is of order \G/K\ > \G/S\ — 4. The core K is the kernel of the 
natural homomorphism </): G — >■ 6*4 which describes the permutation action of G on G/S, and so 
G/K = img(/) < 5*4 is a transitive subgroup of order greater than 4. The only possibilities are A4, 
Z?8 and 54. Because we assume that S is self-normalising, the case G/K = D% is not feasible. The 
statement follows from /3(^4)/|A4| — /3(S'4)/|S'4| = 3/2 (see table [iJ like in the proof of lemma 
10 □ 
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5. Triple Product Property Test Algorithms 

At the moment there are five algorithms to test the TPP. In this section we wiU present them and 
compare their running time. 

We start with two algorithms that directly came from the TPP definition. We define L :— 
Q{S)Q{T)Q{U) as a list, not as a set. Then we count, how many I's are in it. The TPP is fulfilled 
iff there is only one 1 in L. 

Algorithm 5.1. TPPTestNaivC S, T, U ) 

OUTPUT: TPP fulfilled: true / false 

L := Q(s)q(T)q(u) 

if ( #{ 1g € L } > 1 ) then 
return false; 

fi; 

return true ; 

The algorithm above is very naive, because we use the complete list L. This needs a lot of time 
and memory. It is a better idea to search element- wise for I's. 

Algorithm 5.2. TPPTest( S, T U ) 

OUTPUT: TPP fulfilled: true / false 

Qs:=q(S), Qt:-Q(T), Qu := Q(U) 

for s G Qs do 

for t G Qt do 

for u G Qu do 

if ( stu = 1 and ( s =/= 1 ortT^l oru^^l ) ) then 
return false; 

fi; 

od; od; od; 

return true ; 

From Hendrik Orem we know the following equivalent form of the TPP: 

Lemma 5.3. jOrem 20091 Thm. 2.1] Subsets S, T and U of G satisfy the TPP iff 

\S-^\ ■ \U\ = \S-'^U\ and {S-\Q{T)\1)U) D S-^U ^ d). 

Note that for a TPP test based on the lemma above we only need to compute one of the right 
quotients, instead of all three right quotients in the original TPP definition. 

Algorithm 5.4. TPPTestOremC S, T, U ) 

OUTPUT: TPP fulfilled: true / false 

Si '.— S , ^±u '~ SiU 

if ( |Si| • |U| = |Siu| ) then 

if( Si(Q(T) \ l)un Siu = ) then 
return true ; 
fi; fi; 

return false; 

The fourth algorithm comes from the TPP reformulation of theorem |3.1[ it requires a basic 
triple (1 G S* n T n J7) as input: 

Algorithm 5.5. TPPTestMurthy ( S, T, U ) 

OUTPUT: TPP fulfilled: true / false 

Qt Q(T), Qu Q(U) 

if ( Qt n Qu = 1 ) then 

Qs := Q(S) 

if ( Qs n QtQu = 1 ) then 
return true ; 
fi; fi; 

return false; 
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Now we focus on TPP tests for subgroups. We start with the TPP test inspired by theorem 
13.31 about subtransversals. 

Algorithm 5.6. TPPTestMurthyHedtkeGRP( S, T, U ) 

OUTPUT: TPP fulfilled: true / false 
if( SnT = l and SnU=l ) then 
for Xe (S\G)\{S} do 

if ( |XnT| + |Xnu| > 1 ) then 
return false; 
fi; od; 
else 

return false; 
fi; 

return true ; 

In the case where S, T and U are subgroups, we have Q{S) — S, Q{T) = T and Q{U) — U. 
Therefore the remaining test algorithms for subgroups are: 

Algorithm 5.7. TPPTestNaivGRP( S, T, U ) 
OUTPUT: TPP fulfilled: true / false 

L:=S-T-U 

if ( #{ 1g e L } > 1 ) then 

return false; 
fi; 

return true ; 

Algorithm 5.8. TPPTestGRP( S, T, U ) 
OUTPUT: TPP fulfilled: true / false 

for s G S do 

for t e T do 

for u G U do 

if ( stu = 1 and ( s =/= 1 ort^^l oru^^l ) ) then 

return false; 
fi; 

od; od; od; 

return true ; 

Algorithm 5.9. TPPTestDremGRP( S, T, U ) 

OUTPUT: TPP fulfilled: true / false 

Su := S • U 

if ( |S| • |U| = |Su| ) then 

if ( S(T\ l)un Su = ) then 
return true ; 

fi; fi; 

return false; 

Algorithm 5.10. TPPTestMurthyGRP( S, T, U ) 
OUTPUT: TPP fulfilled: true / false 
if ( TnU= 1 ) then 
if ( S n T • U = 1 ) then 
return true ; 

fi; fi; 

return false; 

Now we compare the running time of the different TPP tests. For this we use the search algo- 
rithms described in the next section with the different test algorithms. We use a GAP ( [GAP 2008) ) 
implementation of the algorithms in this paper with the SONATA package ( [Aichinger et. al 2008| ) 
for the Subgroups routine. The interested reader can get the GAP codes from the first author via 
email. 
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Figure 1. Cumulative running time relative to TPPTestNaivGRP of the brute- 
force search for subgroups with the five different test algorithms. (The vertical 
axis is measuring the cumulative runtime of a sequence of 10 runs for each of the 
five different TPP verification algorithms, averaged by 10, so the closer the lines 
are to the zero line, the better.) 



For the subgroup test algorithms we run the brute-force search algorithm 10 times (except for 
the groups of order 64 and 96, because there are more than 200 such groups) for each nonabelian 
group of order less than 128. The results in figure [T] show the mean cumulative running time 
of the search algorithm for all nonabelian groups up to a given order. We recommend to use 
TPPTestMurthyGRP. 

For the subset test algorithms we do the same, but only up to group order 20 and only with 2 
repeats to build the mean value. The results are shown in figure [2] (left). It is very obvious, that 
compared to the naive TPP test the algorithms |5.2[ [5^ and [53] are very fast. Therefore figure [2] 
(right) shows the details. Again, we recommend to use TPPTestMurthy. 




Figure 2. Left: Cumulative running time relative to TPPTestNaiv of the brute- 
force search for subsets with the four different test algorithms. Right: Cumulative 
running time of TPPTestOrem and TPPTestMurthy relative to TPPTest. 



SEARCH AND TEST ALGORITHMS FOR TPP TRIPLES 



9 



6. Brute-force Search 

In this section we describe the brute-force search algorithms that we use for our computations. 
We use the observations from section [2] to reduce the search space for S,T,U C G and the TPP 
tests from the last section. Although it is relatively quick to search for subgroup TPP triples, they 
do not fully describe the TPP capacity of a group, as the following theorem shows. 

Theorem 6.1. For every finite group G, /3{G) > /3g(G) holds. There are groups with j3{G) > 
/3g(G). 

Proof. The first statement is trivial, because the search space for /3 includes the one for /3g. For 
the second statement consider the group Dio = (d, s : = = 1, sds = d~^) of order 10. From 
tablejljwe know that Pg{Dio) = 10. But the subsets S (s), T := {d, s} and U := {1, sd, d^} of 
Dio realize the problem (2, 2, 3) of size 12, and so (3{Dio) = 12. □ 

From the result above we see, that it is necessary so search for subset TPP triples. Therefore 
we present search algorithms for either subgroup TPP triples and for subset TPP triples. 

6.1. Subgroup Search. We use the following algorithm to search for subgroup TPP triples in 
a given nonabelian group. Note that we use the Subgroups routine from the SONATA package 
and that the subgroups returned from that command are ordered ascending by the sizes of the 
subgroups. Our command NonnormalSubgroups is only a filter for the nonnormal Subgroups of 
G. 

Algorithm 6.2. TestGRP( G ) 



OUTPUT: subgroups 



that realize a problem of size /i^g(G) 



(observation 
(observation 
(observation 



3.6) 



2.7) 



2.1) 



(observation 2.5) 



(observation 2.10) 



(observation 2.5) 



1 Umax • ~ I G I , Snax-~G, T^ax ■ ~"[ Iq^" > Umax-~'ClG}' 

2 C : =NonnormalSubgroups (G) ; a:=|C|; 

3 lastS:=| { cGC : |c| < £(G) } I + 1; 

4 for i in [a, a-1, lastS] do 

5 S:=C[i]; 

6 for j in [i-1, i-2 1] do 

7 T:=C[j]; 

8 if ( S n T = 1 ) then 

9 for k in [j-1, j-2, 1] do 

10 U:=C[k]; 

11 if ( |S|-(|T| + lUl - 1) < IGI ) then 

12 if ( ISI-ITI-IUI > n„,ax ) then 

13 if ( T n U = 1 and S n U = 1) then 

14 if TPPfulf illed(S,T,U) then 

15 nmax: = lSl • |T| • |U| ; Smax:=S; Tmax:=T; U^ 

16 break; 

17 fi; fi; 

18 else 

19 break; 

20 fi; fi; od; fi; od; od; return S^^x , T^^t , Umax; 

For the command TPPfulf illed in line 14 of the code we use one of the TPP tests for subgroups 
described in the previous section. The break command in line 16 is used because all other U that 
occur in the inner for loop have a size smaller or equal than J/max and we don't need to test them, 
because we will not get a bigger TPP triple. The break command in line 19 is used because the 
multiplicative size of all other TPP triples with an U from the inner loop have a size smaller than 



(observation 2.12) 
(observation 2.6) 
(observation 2.10) 



=U; 



6.2. Subset Search. The search for subset TPP triples is very similar to the subgroup search. 
But instead of the Subgroups routine we use a method to generate subsets of G. We followed an 
idea of E. Burnett (see [Burnett 2010| ) to generate the subsets on the fly in our code. (We also 
had the idea to generate a binary representation of the subsets. Each subset of G can be identified 



10 



IVO HEDTKE AND SANDEEP MURTHY 



as an element of {0, l}!*^!. But it needs too much memory to generate the whole object {0, 
and therefore we are glad, that we found a way to generate the subsets one after another.) 

Algorithm 6.3. GAP implementation of the GenerateSubset method from E. Burnett 
GenerateSubset := function ( OrderG , SizeS , NumberS ) 

local BoolVector, Offset, NumberZeros , Low, High; 
BoolVector:=[] ; Dffset:=0; NumberS : =NumberS-l ; 
while true do 

NumberZeros : =0 ; Low:=0; High :=Binomial (OrderG- 1, SizeS-1) ; 

while (Numbers >= High) do 

NumberZeros : =NumberZeros+l ; Low : =High ; 

High:=High+Binomial(OrderG-NumberZeros-l, SizeS-1) ; 
od; 



Add (BoolVector , Of f set+NumberZeros+1) ; 



if (Sizes = 1) then 
return BoolVector; 
else 

OrderG : =OrderG-NumberZeros-l ; SizeS : =SizeS-l ; NumberS : =NumberS-Low ; 
Offset : =0f f set+NumberZeros+1 ; 

fi; 

od; 

end; 

The differences to the subgroup search algorithm are: First we compute the set W of all possible 
values of \S\, \T\ and \U\, that means W :— {2,...,\G\ — 1}^. Then we use all observations 
from section [2] to reduce the search space. Now we order the elements of W decreasing by their 
multiplicative size (we use this to stop the algorithm in the case that we found a TPP triple, 
because it is the biggest TPP triple of G). In the outer for loop we go through all i £ W and 
generate in the inner loops the ('^bi^) subsets S that contain the identity of G and the U's and 
T's in the same way. 
for i in W do 

for NrS in [1 . .Binomial (Size (G)-l , i [1] -1)] do 
IndexS := GenerateSubset (Size (G) , i [1] , NrS) ; 
for NrT in [1 . .Binomial(Size(G) -1 , i [2] -1)] do 
IndexT : = GenerateSubset (Size (G) , i [2] , NrT) ; 
if Size(Intersection(IndexS, IndexT)) = 1 then 
for NrU in [1 . .Binomial (Size (G)-l , i [3] -1)] do 
IndexU := GenerateSubset (Size (G) , i [3] , NrU) ; 
if (Size (IntersectiondndexS , IndexU) ) = 1 

and Size(Intersection(IndexT, IndexU)) = 1) then 
S := AsList(G) {IndexS}; ... 

Of course one has to use a TPP test for subsets (see the previous section) . 



7. Selected Results 



In this section we present some computational results of the brute-force search with subsets and 
subgroups. The values in the following tables are: 

• GAP Id: The IdSmallGroup in the SmallGroups Library. 

• £>3(G) resp. D3: The 3-character capacity (see definition 1.4). 

• (3{G) resp. /?: The TPP capacity (see definition 1.3 1. 
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HAP TH 


Group 








P/ J-'3 






Pi fn) 


[6, 1] 


^3 


10 


8 


8 


0.8 


1 00000 

1.33333 


2, 2 


2 


[8,3] 


Ds 




12 


8 


8 


0.666667 


1 


2,2 


2 


[8,4] 


Qs 




12 


8 


8 


0.666667 


1 


2, 2 


2 


[10, 1] 


Dio 


18 


12 


10 


0.666667 


1.2 


3,2 


2 


[12, 1] 


Cs X C4 




20 


16 


12 


0.8 


1.33333 


4,2 


2 


[12,3] 


Ai 




30 


18 


18 


0.6 


1.5 


3,3 


2 


[12,4] 


D12 




20 


16 


16 


0.8 


1.33333 


4,2 


2 


[14, 1] 


Dm 


26 


16 


14 


0.615385 


1.14286 


4,2 


2 


[16,3] 


{Ci X C2 


) X C2 


24 


16 


16 


0.666667 


1 


4,2 


2 


[16,4] 


C4 X C4 




24 


16 


16 


0.666667 


1 


4,2 


2 


[16,6] 


X C2 




24 


16 


16 


0.666667 


1 


4,2 


2 


[16, 7] 


Die 




28 


20 


16 


0.714286 


1.25 


5,2 


2 


[16,8] 


QDlQ 




28 


16 


16 


0.571429 


1 


4,2 


2 


[16,9] 


Q16 




28 


16 


16 


0.571429 


1 


4,2 


2 


[16, 11] 


C2 X Ds 




24 


16 


16 


0.666667 


1 


4,2 


2 


[16, 12] 


C2 X Qs 




24 


16 


16 


0.666667 


1 


4, 2 


2 


[16, 13] 


(C4 X C2 


) X C2 


24 


16 


16 


0.666667 


1 


4, 2 


2 


[18, 1] 


D18 




34 


24 


24 


0.705882 


1.33333 


6, 2 


2 


[18,3] 


C3 X 53 




30 


24 


24 


0.8 


1.33333 


6,2 


2 


[18,4] 


C| X C2 




34 


24 


24 


0.705882 


1.33333 


6,2 


2 


[20, 1] 


C5 X C4 




36 


24 


20 


0.666667 


1.2 


6,2 


2 


[20, 3] 


C5 X C4 




68 


32 


32 


0.470588 


1.6 


4,4 


2 


[20, 4] 


D20 




36 


24 


20 


0.666667 


1.2 


6,2 


2 


[21, 1] 


C7 X C3 


57 


27 


27 


0.473684 


1.28571 


3,3 


3 


[22, 1] 


D22 


42 


28 


22 


0.666667 


1.27273 


7,2 


2 


[24, 1] 


C3 X Cs 




40 


32 


24 


0.8 


1.33333 


4,4 


2 


[24, 3] 


SL2 F3 




54 


36 


36 


0.666667 


1.5 


4,3 


3 


[24, 4] 


C3 X Qs 




44 


32 


24 


0.727273 


1.33333 


4,4 


2 


[24, 5] 


Ci X 53 




40 


32 


32 


0.8 


1.33333 


4,4 


2 


[24, 6] 


D2i 




44 


32 


32 


0.727273 


1.33333 


4,4 


2 


[24, 7] 


C2 X (C3 


X C4) 


40 


32 


24 


0.8 


1.33333 


4,4 


2 


[24, 8] 


( Ce X C2 


) X C2 


44 


32 


32 


0.727273 


1.33333 


4,4 


2 


[24, 10] 


C3 X 7^8 




36 


24 


24 


0.666667 


1 


6,2 


2 


[24, 11] 


C3 X Qs 




36 


24 


24 


0.666667 


1 


6,2 


2 


[24, 12] 


Si 




64 


36 


36 


0.5625 


1.5 


4,3 


3 


[24, 13] 


C2 X Ai 




60 


36 


36 


0.6 


1.5 


4,3 


3 


[24, 14] 


C2 X C2 


< S3 


40 


32 


32 


0.8 


1.33333 


4,4 


2 



Table 1. Computational results for all nonabelian groups of order less than 25. 
(If /? and /3g differ, the corresponding values are printed bold.) 



• /3g(G) resp. (3^: The TPP subgroup capacity (see definition 1.3). 

• /3/lG] resp. /3g/]G]: The so called TPP ratio or TPP subgroup ratio, resp.: A parameter 
that measures the size of the realized problem in relation to the group size. 

• {n,p,m): The parameters of the biggest realized problem. 

7.1. Computational Results for some Small Groups. Our search algorithms allow us to 
compute all the relevant information about TPP triples of subsets for all nonabelian groups of 
order less than 25. The results are shown in table [T] (computing TPP triples of subsets in groups 
of order greater than 24 is too time consuming at present). 

7.2. Computational Results for some 2-Groups. An interesting point is, that the only groups 
we found with the brute-force search for subgroups that achieve /3/-D3 > 0.8 are 2-groups. The 
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HAP Trl 


Group 








Pg/|Cx| 


(n,p, m) 


[64, 226] 


^8 


144 


128 


0.888889 


2 


8,4,4 


[128,29] 


(Ca X (Cg X C2)) X Ci 


304 


256 


0.842105 


2 


16,4,4 


[128,1135] 


(C| X Dg) X C2 


304 


256 


0.842105 


2 


16,4,4 


[128,1142] 


(C| X Dg) X C2 


304 


256 


0.842105 


2 


8,8,4 


[128,1165] 


(C| X Dg) X C2 


304 


256 


0.842105 


2 


8,8,4 


[128,2194] 


C2xDi 


288 


256 


0.888889 


2 


16,4,4 


[128,2213] 


(C2 X C4 X Ds) X C2 


288 


256 


0.888889 


2 


16,4,4 



Table 2. Computational results for selected 2-groups. 



results for groups of order 64 and 128 are shown in table[2] (The results for all nonabelian groups of 
order 256 are still in computation.) Details about the subgroup search can be found in subsection 

EH 

7.3. Selected Results for SL„ and PSL2 F^. We also tested some groups of type SL„ F^ and 
PSL2Fg. The results are shown in the tables [3] and [4[ These are groups that realize a relatively 
high TPP subgroup ratio p{G) :— (i{G)/\G\. No other tested groups obtain ratios bigger than 4. 
Furthermore we present a theoretical result about a lower bound for the TPP capacity of SL3 F2: 

Lemma 7.1. The finite special linear group SL3 F2 realizes (8, 7, 7) via a TPP triple of subgroups, 
and its TPP capacity has the lower bound /3(SL3F2) > 392 = || SL3F2I. 

Proof. We write G := SL3 F2 and define a proper subset S" C G by 



S := 



I X y 
1 z 
1 



x,y,z e F2 



This is a subgroup of order 8, consisting of all the upper unitriangular matrices in G. Define two 
proper cyclic subgroups T,U < G hy 



T:= to := 



111 
1 
1 1 



, U := (uo:= 



1 

1 1 
1 1 1 



The generators to and uq, as defined above, have order 7 in G, so T and U are cyclic of order 7. 
We know that T O U = 1 (otherwise to would generate U, and we would have T = U, which is 
not true). In TU :~ {^0^0 ■ ^ ^7} C G all nontrivial products IqUq, with i,j G Z7 not both 
zero, are non- upper triangular matrices, and therefore T, U,TU C (G\ S*) U 1. Therefore, we have 
SnrU ^1. Thus (S, T, U) is a (basic) TPP triple of subgroups of G, and by it G reahzes (8, 7, 7) 
of size 8 • 7 • 7 = 392. The TPP capacity /3(G) is at least its TPP capacity via subgroups /3g(G), 
which, by this result, is at least 392. Since SL3F2 = PSL3F2 = PSL2F7, this result is also true 
for these latter groups. □ 

Note that it follows from the tables |3] and [I] that /3g(SL3F2) = 392 holds. 

7.4. Further Results. At the moment many computations are still running on three supercom- 
puters at the universities of Jena and Halle- Wittenberg. A first goal is to compute the TPP 
subgroup capacity /3g for all nonabelian groups of order up to 1000 (except 512 and 768). We 
also hope that there exists an algorithm that can effectively produce basic subsets, because the 
computation of the TPP subset capacity /? has a complexity of O (81*^1) without such an algorithm 
and it takes 8 hours to compute f3 for one group of order 24. All results (that means of the subset 
search for groups of order up to 24 and of the subgroup search of groups up to order 1000) can be 
found at http://www2.informatik.uni-lialle.de/da/hedtke/tpp/. Note that some results are 
still missing, because they are in computation at the moment. 
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Group 




|G| 


/3g(G) 




^g/|G| 


(n, p, m) 


PSL2 F2 


10 


6 


8 


0.8 


1.33333 


2,2,2 


PSL2F3 


30 


12 


18 


0.6 


1.5 


3,3,2 


PSL2F4 


244 


60 


108 


0.442623 


1.8 


12,3,3 


PSL2 F5 


244 


60 


108 


0.442623 


1.8 


12,3,3 


PSL2 F7 


1126 


168 


392 


0.348135 


2.33333 


8,7,7 


PSL2F8 


4072 


504 


1372 


0.336935 


2.72222 


14, 14, 7 


PSL2 Fg 


3004 


360 


972 


0.323569 


2.7 


12,9,9 


PSL2Fn 


7038 


660 


1980 


0.28133 


3 


55,6,6 


PSL2 F13 


13556 


1092 


3276 


0.241664 


3 


78, 7, 6 


PSL2F17 


40252 


2448 


10368 


0.257577 


4.23529 


24,24,18 


PSL2F19 


63646 


3420 


14400 


0.226251 


4.21053 


60,20,12 



Table 3. Computational results for selected groups of type PSL2Fg. 



Group 




\G\ 


/3g(G) 






(n,p, m) 


SL2 F2 


10 


6 


8 


0.8 


1.33333 


2,2,2 


SL2 F3 


54 


24 


36 


0.666667 


1.5 


4,3,3 


SL2 F4 


244 


60 


108 


0.442623 


1.8 


12,3,3 


SL2F5 


540 


120 


216 


0.4 


1.8 


24,3,3 


SL2 F7 


2198 


336 


784 


0.356688 


2.33333 


16,7,7 


SL2 Fs 


4072 


504 


1372 


0.336935 


2.72222 


14,14,7 


SL3F2 


1126 


168 


392 


0.348135 


2.33333 


8,7,7 



Table 4. Computational results for selected groups of type SL„Fg. 
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